TUTP-96-2 
July 1, 1996 
Revised Oct. 2, 1996 



RESTRICTIONS 
ON NEGATIVE ENERGY 
DENSITY 
IN FLAT SPACETIME 



L.H. FordF] and Thomas A. RomanQ 

Institute of Cosmology 
Department of Physics and Astronomy 
Tufts University 
Medford, Massachusetts 02155 

Abstract 

In a previous paper, a bound on the negative energy density seen by an 
arbitrary inertial observer was derived for the free massless, quantized scalar 
field in four-dimensional Minkowski spacetime. This constraint has the form 
of an uncertainty principle-type limitation on the magnitude and duration of 
the negative energy density. That result was obtained after a somewhat com- 
plicated analysis. The goal of the current paper is to present a much simpler 
method for obtaining such constraints. Similar "quantum inequality" bounds 
on negative energy density are derived for the electromagnetic field, and for the 
massive scalar field in both two and four-dimensional Minkowski spacetime. 
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1 Introduction 



In quantum field theory, unlike in classical physics, the energy density may be un- 
boundedly negative at a spacetime point. Such situations entail violations of all 
the known classical pointwise energy conditions, such as the weak energy condition 
This fact has been known for quite sometime 0. Specific examples include the 
Casimir effect [|3], Q and squeezed states of light ||, both of which have observational 
support. The theoretical prediction of black hole evaporation || also involves neg- 
ative energy densities and fluxes in a crucial way. On the other hand, if the laws 
of quantum field theory place no restrictions on negative energy, then it might be 
possible to produce gross macroscopic effects such as: violation of the second law of 
thermodynamics 0, |8| or of cosmic censorship & 10 , traversable wormholes [11] , [L2| , 
'warp drive" 0], and possibly time machines |J, 14|. As a result, much effort has 



been recently directed toward determining what constraints, if any, the laws of quan- 
tum field theory place on negative energy density. One approach involves so-called 
'averaged energy conditions" (see, for example, Jl5|]-[l9[j), i.e., averaging the local en- 



ergy conditions over timelike or null geodesies. Another method employs "quantum 
inequalities" (QI's) JF] p0| , which are constraints on the magnitude and duration of 
negative energy fluxes and densities. The current paper is another in a series which 
is exploring the ramifications of this approach [PT|-||25||. (For a more comprehensive 



discussion of the history of these topics, see the introductions of Refs. |22| and the 
references therein.) 

The QI's have the general form of an inverse relation between an integral involving 
the the energy density or flux over a finite time interval and a power of that interval. 
More precise forms of the inequality were originally derived for negative energy fluxes 
|2U|| , and later for negative energy density |^T], B^|. This form of QI's involves "folding" 
the stress energy tensor into a "sampling function", i.e., a peaked function of time 
whose time integral is unity. For example, it was shown in Ref.[[T[| that for the free 
quantized massless scalar field in four-dimensional Minkowski spacetime, 



(T 00 )dt 



vr 7-oo t 2 + t 2 ~ 32vr 2 t n 4 



for all choices of the sampling time, t^. Here (Too) is the renormalized expectation 
value of the energy density evaluated in an arbitrary quantum state in the frame 
of an arbitrary inertial observer whose proper time coordinate is to. The physical 
implication of this QI is that such an observer cannot see unboundedly large negative 
energy densities which persist for arbitrarily long periods of time. The QI constraints 
can be considered to be midway between the local energy conditions, which are ap- 
plied at a single spacetime point, and the averaged energy conditions which are global, 
in the sense that they involve averaging over complete or half-complete geodesies. The 
QI's place bounds on the magnitude and duration of the negative energy density in 
a finite neighborhood of a spacetime point along an observer's worldline. 

These inequalities were derived for Minkowski spacetime, in the absence of bound- 
aries. However, we recently argued [23] that if one is willing to restrict the choice 
of sampling time, then the bound should also hold in curved spacetime and/or one 
with boundaries . For example, we proved that the inequality Eq. (jl]) holds in the 
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case of the Casimir effect for sampling times much smaller than the distance between 
the plates. It turns out that this observation has some interesting implications for 
traversable wormholes [E3]. Quantum inequalities in particular curved spacetimes, 
which reduce to Eq. ([I]) in the short sampling time limit, are given in Ref. |[24|| . 

In the original derivation of Eq. ([I]), we used a rather cumbersome expansion of 
the mode functions of the quantum field in terms of spherical waves. The goal of the 
present paper is to present a much more transparent derivation of QI bounds, based on 
a plane wave mode expansion. In so doing, we prove new QI constraints on negative 
energy density for the quantized electromagnetic and massive scalar fields. In Sec. 0, 
we derive a QI bound for the massive scalar field in both four and two-dimensional 
Minkowski spacetime. Our earlier result, Eq. ([[]), is recovered as a special case when 
the mass m goes to zero. A similar bound is obtained for the electromagnetic field in 
Sec. H Our results, and their implications for the existence of traversable wormholes, 
are discussed in Sec. f|. Our metric sign convention is (— , +, +, +). 



2 The Massive Scalar Field 
2.1 Four-Dimensions 

In this section we derive a Ql-bound on the energy density of a quantized uncharged 
massive scalar field in four-dimensional flat spacetime. The wave equation for the 
field is 

(□-m 2 )0 = O, (2) 

where □ = , rf v d^d v . We can expand the field operator in terms of creation and 
annihilation operators as 

= EK/k + 4/ k *). (3) 

k 

Here the mode functions are taken to be 
where 



lu = \/|k| 2 + m 2 , (5) 

m is the rest mass, and V is the normalization volume. The stress tensor for the 
massive scalar field is 

T,u, = - + m 2 (f) 2 ) . (6) 

The renormalized expectation value of the energy density, in an arbitrary quantum 
state is 

IT \ - Re \" + k ' • k ) [/+ v _i( w '-w)t , /. _ \_-i(w'+w)t 



'UO'UJ 
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Here the energy density is evaluated in the reference frame of an inertial observer, 
at an arbitrary spatial point which we choose to be x = 0. The time coordinate t is 
the proper time of this observer. Because of the underlying Lorentz invariance of the 
field theory, we may of course choose the frame of any inertial observer in Minkowski 
spacetime. 

Following Refs. [^0]. |2l| , we multiply (Too) by a peaked function of time whose 
time integral is unity and whose characteristic width is to- A convenient choice of 
such a function is t /[7i(t 2 + Define the integrated energy density to be 



P 



tn 



71 



(Tqq) dt 
t 2 + tl • 



(8) 



Substitution of Eq. (0) into Eq. (H) and performance of the integration yields 
P 



Re ^ (u'u + k' • k) r 



Re ^ 

+ — E- 



m 



UJ'UO 



(aUk) e-^'-^ t0 + (a k /a k ) e'^^ 
(4,a k ) e ->'-»l*> - (a k ,a k ) e -(<*'+<*)<° 



(9) 



Now write k' ■ k = k' x k x + k' y k y + k' z k z , and apply the lemma in Appendix B to each 
of the sums in Eq. (0), e.g., 



k',k 

We then obtain 



-= «,a k ) e 



W-u\to > ( a t, ak ) e -(-'+^ 

k',kVW'w 



(10) 



> Re ^ (a/a; + k' ■ k) (w , +ul)to 



(a k ,a k ) + (a k /a k ) 



+ 



Re . , m 2 
2^kVkV^ 



-(u/+u;)to 



(a k ,a k ) - (a k /a k ) 



'11^ 



For the first sum (i.e., the one involving u'u), apply the first lemma in Appendix 
A, where we take hj = (l/2)(Ju /V) e~ w<0 . Apply the same lemma to each term in 



the sum which involves k' • k, taking hj = (ki/2y/uV) e wt °, with i = x,y,z. For 
the last sum, involving m 2 , use the second lemma in Appendix A choosing hj = 

{m/2VujV)e- ultQ . The result is 



> 



2V 



E 



(12) 



where we have used Ikl 
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m 2 . Now let V — > oo, in which limit TJ k 



(V^/87r 3 ) /f^ d 3 k. Performing the angular integrations and changing variables we get 

1 



P> 



m 2 uj 2 e~ 2ujta 



(13) 



G(y) 




2 4 6 8 10 

y 

Figure 1: 
The graph of G(y) versus y. 



Let x = 2uto, y = 2mt - We may then rewrite the integral in Eq. (|13|) as 



P>~ aA 2+4 dx^/x 2 -y 2 x 2 e~ x . (14) 

D47T tg Jy 



If we now let 



1 C°° I 

G{y) = -J dx^/x 2 -y 2 x 2 e- x , (15) 
then our Ql-bound may be written as 

^-32^°^' <16) 

for all to- A plot of G(y) versus y is given in Fig. 1. 

We see that asm^O (for fixed t ), y — > and G(y) — > 1, hence in this case our 
bound reduces to 

^-32^V' (17) 

for all to, which is the QI for the massless scalar field - originally obtained in [^TJ by 
a much more complicated method. As m — > oo (for fixed t ), y — > oo, and we see 
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from the graph that G(y) — » 0. Note that our bound becomes more stringent as m 
increases, and hence it becomes increasingly difficult to produce large negative energy 
densities. This result is not surprising, since now one has to overcome the positive 
rest mass energy of the field quanta. For i ^ l/ m (with fixed m), corresponding 
to sampling times much larger than the Compton wavelength of the particle, we also 
have y — ► oo and G(y) — ► 0. Note that, due to the factor of G(y), the right-hand 
side of Eq. ([16]) vanishes more rapidly than the right-hand side of Eq. ( |I?D in the 
t — > oo limit. This can be interpreted as showing that negative energy due to a 
massive scalar field must be more highly localized than that due to a massless field. 



2.2 Two-Dimensions 

The derivation of a QI bound on the energy density of a massive scalar field in two- 
dimensional spacetime is very similar to the discussion given in the previous section 
for the four-dimensional case. One follows the same steps as in that treatment, but 
replacing the normalization volume V with a periodicity length L, and k and x by k 
and x, respectively. This procedure yields 

p > ~ E ^~ 2ujt0 , (is) 

k 

where we have used k 2 = uj 2 — m 2 . It is at this stage that the 2D analysis differs from 
the 4D case. Now let L — > oo, in which limit J2k ~> (L/2ir) J^dk. Rewriting the 
integral in terms of u, one obtains 

P>-— / du. (19) 



Let 



2n Jm y/u 2 — m 



oo ^ g— 2^o 



This expression can be rewritten [2£| in terms of modified Bessel functions of the 
second kind, K v {y), as follows 

I = -m 2 K[(2mt ) , (21) 

where the prime denotes differentiation with respect to the argument. We now employ 
the relation pTj 

K[(y) = -\[K {y) + K 2 {y)\ , (22) 
where again y = 2mto- The integral / can then be written as 

1 



y 2 [K (y) + K 2 (y)} 



2 

it 2 ' 

1 F(y). (23) 



At 



2 
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F(y) 




Thus our Ql-bound may be written as 

P>-^F(y), (24) 

for all to- The function F(y) is plotted against y in Fig. 2. 

Asm->0 (for fixed t ), F(y) — > 1, and our inequality becomes 

>*-S5»- (25) 

for all toj which is the QI for the massless scalar field in 2D pifl . For large y, the 
behavior of the graph is qualitatively similar to the 4D case. However, for y in the 
approximate range < y < 1.2, there appears to be a small peak in the value of 
F(y). In the case of flat spacetime, this seems to be an artifact of 2D, since no 
corresponding peak occurs in the 4D case |^?J . Rather than indicating that one could 
actually supercede the m = bound in the 2D massive case for certain values of y, 
this result may perhaps imply that the bound given in Eq. (|2"4D is not optimal over 
this range. In either case, the relative height of the peak above the value of F(y) = 1 
(corresponding to the massless case) is too small to produce a dramatic change in 
the QI bound. As we argued in Ref . |23fl , the constant on the right-hand side of the 
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inequality would typically have to be larger by many orders of magnitude to result 
in large macroscopic effects. 



3 The Electromagnetic Field 

In this section, we derive a Ql-bound for the energy density of the free quantized 
electromagnetic field in 4D Minkowski spacetime. The absence of source charges and 
the imposition of the Coulomb gauge condition imply 

A = $ = , 

V-A = 0, (26) 

where A is the vector potential, and $ is the scalar potential. The wave equation is 

□A = . (27) 

The vector potential can be expanded in terms of creation and annihilation operators 
as 

A ( x > *) = XX a kAekA/k + al x e kX fk) , (28) 

kA 

where the e kA are unit linear polarization vectors and thus A = 1,2. Here the sum 
SkA = Z)k J2x- The mode functions / k are again given by 

h = — ^e i(kx -"' } , (29) 

i.e., we assume periodic boundary conditions in a box with normalization volume V . 
We will adopt the convention that one of the polarization vectors changes direction 
as k — > — k, but the other does not. That is, we have 

ekixe k2 = k, 
e_ki = — e k i > 

e_k2 = e k2 , (30) 

so that e_ki x e-k2 = — k. The commutation relations for the creation and annihila- 
tion operators are 

[a kA , a k , v ] = 5 kk / 5\\> , (31) 
KcA, Ctk'A'] = , (32) 

[4a, a l\>] = • (33) 
With the gauge choice made in Eq. (p6|), the electric field operator becomes 

E = -A = J2iu (a kA e kA / k - <4 A e kA / k *) . (34) 

kA 



8 



Similarly the magnetic field operator is given by 

B = VxA = ^iw (a kX b kX f k - al x b kX fk) 



Here we define 
so that 

and therefore 
and 



-kl 



b kX = kxe kA , 

6 kl = kxe kl = e k2 , 
6 k2 = kxe k2 = -e k i 

-kxe_ k i = e k2 = b 



ki 



b-k2 = -kxe_ k2 = e ki = -6 k2 . 
The stress energy tensor for the electromagnetic field is 

T 00 = ^(E 2 + B 2 ). 



(35) 

(36) 

(37) 
(38) 

(39) 

(40) 

(41) 



If we compute the renormalized expectation value (T o), in an arbitrary quantum 
state \ip), evaluated at x = 0, and fold this into our sampling function we obtain 

. _ t f°° (T 00 )dt 



7T J—oo t' 2 ~\~ tg 

Re 



— PkA • e k 'A') + (b kX ■ 6k' A')] 

^ ' lr 1*-' 



k.k' 

A, A' 



X 



(ai x a k ,y) e-l— + (a kA a k , A ,> e -(«^')*o 



(42) 



Now write e kA • e k / A / = (e^ (e k 'A')a; + ( g kx) y (e k ' A ') y + (^\) z (e k ' A ') z > and similarly 
for the dot product involving the b vectors. Expanding Eq. ( f42|) in this fashion, one 
gets six terms. To each of these terms we apply the lemma in Appendix B, e.g., 



£ V^?(e kX ) x (e k ,y) x (at x a k , x ,) e^'^ 

k,k' 

A, A' 

> X; vW(e kA ) x (e k , v ) ;E (a kA a k , v )e-^') to 



(43) 



k.k' 

A, A' 



Recombining the terms, we then obtain 

Re 



k,k' 

A, A' 



X 



(a kA Q k'A') + (flkAOk'A') 



-(w+u/)to 



(44) 
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An analogous separation of the dot products in Eq. (fPj) into components, repeated 
application of the first lemma in Appendix A, and a recombination of terms results 
in 

P > ~ E ™™* , (45) 

ZV kA 



where we have used hj = (1/2) (e kA )i x Ju/V ' e'^ and hj = (1/2) (6 kA ). x yfu/V e'^ , 
with i = x,y, z. Evaluation of the sum over A gives 

P>-^E^" 2 ^. (46) 



P>-^ duuj z e- 2ult0 . (47) 



Now let V — > oo, so 

o > - 

2vr 

An evaluation of the integral gives us our desired result: 

^-TW' (48) 

for all to- Comparison with the Ql-bound for the massless scalar field, Eq. ([H]), shows 
that the right-hand side of the bound in the electromagnetic field case is smaller (i.e., 
more negative) by a factor of 2. This is just what one would expect due to the two 
polarization degrees of freedom. 



4 Conclusions 

In the current paper, we have derived "quantum inequality" (QI) bounds on negative 
energy density for the quantized (uncharged) massive scalar and electromagnetic fields 
in Minkowski spacetime. The bounds take the form of uncertainty principle-type in- 
equalities on the magnitude and duration of the negative energy density seen by an 
arbitrary inertial observer. (Note however that the energy-time uncertainty principle 
is not used as input to derive the QI bounds.) Originally we had derived such a bound 
for the massless scalar field in both two and four-dimensional Minkowski spacetime 



21|j . Our earlier four-dimensional derivation was performed using a rather compli- 
cated expansion of the mode functions in spherical waves. The present treatment is 
considerably simpler, as the analysis is done in a plane wave mode representation. 
We recover our previous results for the massless scalar field as the m = case of the 
massive scalar field. For the massive scalar field, we find that in general it is more 
difficult to obtain sustained negative energy densities than in the massless case. This 
is not surprising, as now one must overcome the positive rest mass energy. In the case 
of the electromagnetic field, the right-hand side of the bound is slightly weaker, i.e., 
by a factor of 2, than in the massless scalar field case. This is also to be expected, 
given the two polarization degrees of freedom in the former case. For all of our QI 
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bounds on energy density, in the infinite sampling time limit we sample over the en- 
tire timelike geodesic of the observer and obtain the averaged weak energy condition 
(AWEC): 

roo 

T^uV dr > , (49) 



where is the unit tangent vector to the geodesic and r is the observer's proper 
time. Hence, in Minkowski spacetime, the AWEC can be derived from the QI bound. 

Note that in Minkowski spacetime, our bounds hold for an arbitrary choice of 
the sampling time, to- Recently we argued that such bounds should also hold in a 
curved spacetime and/or one with boundaries, if the sampling time is restricted to 
be much smaller than the smallest local radius of curvature and/or the distance to 
any boundaries in the spacetime |p3|. An explicit example of the validity of this 
assumption has been given in Ref. jp4[| . There a QI is derived in the static closed 
and open Robertson- Walker universe models. For a choice of sampling time which is 
small compared to the local radius of curvature, the QI reduces to our flat spacetime 
bound. 



The application of our bound to Morris-Thorne wormhole spacetimes p3|| , in the 
small sampling time limit, implied that typically either a wormhole must be only 
slightly larger than the Planck size, or that the negative energy density must be 
concentrated in an extraordinarily thin band near the throat. As an example, in one 
case for a wormhole with a throat radius of 1 m, the negative energy density must be 
concentrated in a band no thicker than about a millionth of a proton radius. Even 
for a wormhole with a throat radius the size of a galaxy, the band of negative energy 
must be no thicker than about 10 proton radii (for this particular type of wormhole) 
||29|1 . That analysis assumed that the stress energy tensor which maintained the 
wormhole consisted of quantized massless scalar fields. More specifically, the (T^) 
which generates the wormhole geometry was assumed to be the expectation value of 
the stress energy tensor operator of the massless scalar field in a suitable quantum 
state. The results of the present paper indicate that the same restrictions on the 
wormhole length scales will also apply if one tries to make the stress energy of the 
wormhole out of electromagnetic or massive scalar fields. One might hope that the 
constraints imposed by these bounds could be circumvented by the superposition of 
many fields, each of which individually satisfies a QI bound. However, we showed 
|23H that in practice to achieve significant macroscopic effects, one needs either on the 
order of 10 62 fields or a few fields for which the numerical constant on the right-hand 
side of the QI bound is many orders of magnitude larger than in the cases examined 
to date. Neither of these possibilities seem very likely. It therefore appears probable 
that Nature will always prevent us from producing gross macroscopic effects with 
negative energy. 
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Appendix A 

In this appendix we will establish two lemmas on sums of the expectation values of 
products of creation and annihilation operators. The discussion is a slight generaliza- 
tion of the argument presented in the Appendix of Ref . |TD| . The idea for this method 



of proof was originally suggested to us by Flanagan pi . The first lemma was proven 



by a more complicated argument in Appendix A of Ref. 

Lemma 1: 

Define the operator Oj by 

4 = EM% + 4> (Al) 

j 

where j is a generalized mode label, and the h/s are assumed to be real. Then it is 
easily seen that Oj is a Hermitian operator. We may form the expectation value 

(Opj) > . (A2) 

This follows from the fact that the left-hand side is simply the norm of the state 
vector = Oj\ijj), where \ip) is the quantum state in which the expectation value 
is taken. We may expand this expression as 

(Ofii) = E MiK4 a i> + ( Q K> + + (4 a i)l + E ( h if- ( A3 ) 

3,3' 3 

Use the fact that (aj/ a i) = ( a ] a j') to write 

(Ofij) =2ReJ2 + (Oj'Oj)] + E ( h 3-f- (A4) 

3,3' 3 

It follows from Eq. flA"2] ) that 

2 Re E M;[<4%> + ^ "E (hj) 2 . (A5) 

3,3' 3 

Lemma 2: 

Define the operator Pj by 

P^i^hM-aj), (A6) 

3 

where j is a generalized mode label, and the h^s are again assumed to be real. Then 
the operator Pj is also Hermitian. Utilizing the same chain of reasoning used to 
establish Lemma 1, one can show that 

2 Re h j' h A( a ]' a j) ~ (%'%)] > ~E ( h j) 2 - (A7) 

3,3' 3 
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Appendix B 

In this appendix, we wish to prove a lemma which generalizes the inequality proven 
in Appendix B of Ref . PU| . Consider the sum 



S m = $>>i' ( a W) e-l^^'l* . (Bl) 



Here the summation is over a finite set of modes with frequencies Uj, with uj m = 
maxjcuj}, and pj is an arbitrary complex function of the mode label j. (In Ref. [20] , 
it was assumed that the Uj oc j. That assumption is not needed here.) We first note 

that S m is real as a consequence of the fact that (ajaj/)* = (a],aj). Next we define 
the sum 

Sm = Y,P>j' («]%') e-^ + ^ t0 , (B2) 

33' 

which is also real. We wish to prove that 

S m — S m . (B3) 

First note that 

S m > . (B4) 
This follows from the fact that any sum of the form 

^2PjPj'(<4 a f) ( B5 ) 
jj' 

is non-negative as a consequence of its being the norm of the state vector 

l*> = E^#> (B6) 

j 

where is the quantum state in which the expectation values in Eqs. ( |B1| ) and (|B2| ) 
are to be taken. 

Let us assume that the eigenfrequencies are discrete, as will be the case with a 
finite quantization volume, and that the lowest frequency is greater than zero. We 
can then order the modes in increasing frequency so that < u < uj\ < ■ ■ ■ < uj m . 
In general, there can be several modes with the same frequency. Let I label the 
distinct frequencies, and let n(£) label the various modes all having frequency uj£. 
Then £j = Y,eJ2 n (£)- 

Define 

A jf = P * Pr (a]a f )e-^-^\ (B7) 

and 

B 3j' = Pj V 3 > (a]a f ) e^ + ^' )t0 , (B8) 

so that 

S m - S m = Y1( A 33' - B jj') ■ ( B9 ) 

jj' 
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Note that 

ujj + u)ji — \ujj — u)j>\ = 2 Vtjji , (BIO) 
where IX, y is the lesser of Uj and uy. We can then write 

Ajji - B jf = p* Pj > (a]a r ) (e 2 ^' *° - 1) e~^ + ^ to . (Bl 1) 

If we let gj = e _a;j *° pj and define 

C ee = E 9*j 9f («K> , (B12) 

n{t) 
n(t') 

and 

D ee = (e 2n n' t0 -1)C U/ , (B13) 

then we have that 

m m 

S m -S m = - B jf ) = E (e 2n "' t0 - l)C te = E £W . (B14) 

Our goal will be to prove that the right-hand side of this expression is non-negative. 

As a prelude, let us note that any sum of the form J2ei' Cee, where I and £' range 
over the same set of values, is non-negative as a consequence of its being the norm of 
a state vector. In particular, C rnrn > 0. The quantity Y^™e=oDee is the sum of all of 
the elements of an (m + 1) x (m + 1) matrix. Our plan is to prove the positivity of 
this sum by working our way through the matrix, beginning in the lower right-hand 
corner. The element in this corner is 

D mm = (& lUmto -l)C mm >0. (B15) 

Next consider the 2x2 matrix formed by the four elements in this corner. The sum 
over these elements is 

D m — l,m— 1 ~l~ -^m— l,m ^ra,m-l D m m 

— (e 2Wm ~ lto — l)(C m _i )?n _i + C m _i )?n + C mjm _i) + (e 2a,m *° — i) C mm 

^ (6 l)(Cm— l,m— 1 "I - Cm— l,m "I - ^m,m-l "I - Crnrn) 

> 0. (B16) 

Here we have used the facts that C mm > 0, C m _i jm _i + C m - hm + C m , m _i + C mm > 0, 
and that e 2Wmt ° - 1 > e 2Wm - lt() - 1. 

Now suppose that we have established that 

m m 

E > ( e2u;i+lt0 - !) E C«'>0- (B17) 

We wish to show that this implies that 

m m 

E D ee > (e 2oJLt0 - 1) £ CW > . (B18) 
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The additional terms which are added in going from Eq. ( |B17| ) to Eq. ( |B18| ) are those 
which lie in the same row as and to the right of Dll and in the same column as and 
below Dll, as illustrated below: 



£00 
D w 



D, 



01 



Dll 
Dl+i,l 



Dl,l+i 
Dl+i,l+i 



Dltu 
Dl+1,ui 



D 



mL 



D 



m,L+l 



D r 



Note that ue' increases to the right and uj£ increases downward. Thus Qee' = for 
all of these elements. The sum of these terms is 



AD L = D LL + D L , L+ i + ••• + D Lm + D L+hL + V D 



mL 



J2u! L t() 



Now we have that 



,L+1 + \~ C Lm + C L +1,L + • ■ • + C m L) 



(B19) 



ad l + £ D * 



--L+1 



> ADl + (e 2aJL+1 *° — 1) Yl C«' 

e,£'=L+i 



> AD L + (e 2uJLt0 

rr 

= ( e 2 ^°-l) Q 



1) E 

ei'=L+i 



Li'=L 



> 0. 



(B20) 



We have established that Eq. ( pi8| ) holds for L = m and for L — m — 1. We have 
further established that if it holds for one value of L, then it holds for the next smaller 
value of L. It now follows by induction that 



and hence that 



E D w>o 

e,e=o 



[B2D 



(B22) 



Recall that the sums in both S m and S m are over a finite set of modes with max- 
imum frequency uj m . However, we may now take the limit in which both the number 
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of modes and u m become infinite. So long as the resulting sums are convergent, which 
is the case for the systems with which we are concerned, then we have that 

oo oo 

E Pj Pf ( a W) e~^-^' lt0 > £ p* p f (a]a f ) e -^+^ to , (B23) 
jj' jj' 

which is our final result. 
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